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We investigated the spatio-temporal seismicity parameters of shallow (0-70 km) and intermediate (70-300 km) 
depth earthquakes in the Hindu Kush region, which is characterized by the occurrence of large earthquakes in a 
small zone of intense activity. By way of comparison, intermediate depth earthquakes dominate the Hindu Kush 
seismicity. Using a catalogue of 3820 earthquakes, we determined the various earthquake histograms, time- 
series plots, variations of frequency-magnitude distribution (b-value) and seismicity rate changes (z-value). Both 
time periods, encompassing pre- and post- 1964 subsets of earthquakes, differ significantly in terms of reporting 
earthquakes, b- and z-value for shallow and intermediate seismicity. The b-values appear to be lower (1.07 and 
0.90) for intermediate depth earthquakes than do the b-values produced by shallow earthquakes (1.32 and 1.06) 
using the pre- and post-1964 data sets. The three low b-value (< 0.9) areas have been identified within the 
intermediate depth seismicity region of the Hindu Kush. The z-value maps show that the earthquakes with 
magnitude of 7.0 and above can only be seen in the intermediate depth earthquake maps. 


1. Introduction 

The arcuate shaped Himalaya-Karakoram-Hindu Kush-Pamir 

(HKHP) ranges, which collectively comprise the northernmost part of 
Pakistan and adjacent countries, is the most active seismotectonic ter- 
rene (Yin, 2006). These ranges are episodically formed as a result of the 
Hindu Kush, Karakoram, Kohistan Island Arc and the Indian plate col- 
lision with the Eurasian plate following the closure of the Tethys Ocean 
from the Jurassic to Eocene (Faisal et al., 2014). Collectively, the 
northernmost part of Pakistan is the eastern continuation of the greater 
Alpine-Himalayan system, which extends over approximately 7000 km 
from the Mediterranean Sea in the east to the Sumatra arc of Indonesia 
in the west (Kearey and Vine, 1990). This immense orogenic system had 
a profound effect on thrust- and fold-induced topography, tectonic 
burial, global atmospheric circulation, seismicity, geothermal fields 
that control diagenetic/metamorphic and anatectic processes, foreland 
sedimentation, exhumation and structural architecture of the crust that 
influence the socio-economic life of more than one-fifth of the world’s 


population. Seismic activity in Pakistan is concentrated at the con- 
vergent plate boundaries defined by the Rushan-Pshart suture, Wakhan 
suture, Main Karakorum Thrust (MKT) and Main Mantle Thrust (MMT; 
Faisal et al., 2014; Fig. 1). Moreover, the concentration of movements 
along the Main Boundary Thrust, Hazara-Kashmir Syntaxis and Nanga 
Parbat Syntaxis indicate that the interior parts of the Indian plate are 
still tectonically active (Fig. 1; Molnar and Tapponnier, 1975; Kearey 
and Vine, 1990; Treloar et al., 1992; Beck et al., 1995; Kazmi and Jan, 
1997; Burg et al., 2005; DiPietro et al., 2008; Faisal et al., 2014; 
Rehman et al., 2014, 2015). Rehman et al. (2015) divided northern 
Pakistan into three broad seismotectonic domains, known as Himalaya, 
Karakoram and Hindu Kush. In general, shallow earthquakes char- 
acterize the Himalaya and Karakoram domains, while the Hindu Kush is 
famous for intense and intermediate depth earthquakes. The 1974 
Pattan earthquake, at a depth ranging from 8 to 10 km (Jackson and 
Yielding, 1983), and the 2005 Kashmir earthquake, at a depth of 26 km 
(Peiris et al., 2008), are the typical examples of shallow earthquakes in 
the Kohistan Island Arc and Himalayas, respectively. Conversely, the 
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Fig. 1. Map of the present day tectonic grain of Pakistan 
showing plate boundaries, regional thrusts, strike slip faults 
and syntaxial structures (after Kazmi and Jan, 1997; 
Rehman et al., 2015). 
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Fig. 2. Earthquake epicenters (blue and orange) between 1900 and 
2007 from the earthquake catalog of Rehman et al. (2014). Blacklines 
denote faults. Abbreviations of fault names: AN, Andarab; CB, Central 
Badakhsan; CH, Chaman Fault; DZ, Darvaz; GA, Gardez; GF, Ghaza- 
band Fault; HV, Henjvan; MBT, Main Boundary Thrust; MKT, Main 
Karakoram Fault; MMT, Main Mantle; SA, Sarobi; SP, Spinghar 
(modified after Wheeler et al., 2005; NESPAK, 2006). 
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Table 1 

Basic statistics of implemented earthquake catalogues. 


Number 

Catalogue type 

Year 

Reported earthquakes 

Depth (km) 

1 

ISC 

1925-2005 

2178 

1-292 

2 

EHB 

1964-1995 

1141 

1-300 

3 

NEIC 

1994-2007 

229 

7-282 

4 

GSHAP 

1902-1983 

216 

1-300 

5 

SSR 

1962-1989 

42 

1-240 

6 

ABE 

1908-1949 

14 

40-240 


Hindu Kush has mostly experienced intermediate earthquakes. Ex- 
amples of intermediate earthquakes in the Hindu Kush include the 20 
February 1998 earthquake (Mw = 6.4, depth = 236 km), the 3 March 
2002 earthquake (Mw = 7.4, depth = 256 km) and the 26 October 
2015 earthquake (Mw = 7.5, depth = 210 km), which caused damages 
in Afghanistan and Pakistan (USGS, 2002, 2015; Rehman et al., 2015). 
The predominant damages reported by those earthquakes were ca- 
sualties. The 1998 and 2002 earthquakes killed 1 and 113 people, re- 
spectively, while the October 2015 earthquake claimed 280 lives and 
injured 1770 individuals in Pakistan (USGS, 2015). The Malakand re- 
gion of the Khyber-Pakhtunkhwa, located in Pakistan, was affected the 
most in the 2015 earthquake. There were a total of 115 deaths and 538 
injuries reported in Afghanistan in the 2015 earthquake. 


In the Hindu Kush region, the Moho has been reported at depth of 
approximately 70 km and separates the upper crust, which has less 
seismicity, from the relatively high mantle seismic zone (Roecker, 
1981; Ram and Yadav, 1982; Zarifi and Havskov, 2003; Wheeler et al., 
2005). According to Zarifi and Havskov (2003), there is no general and 
unique interpretation for the possible origin of the intermediate 
earthquakes in the Hindu Kush, which could be associated with a 
subduction zone. The concentration of earthquakes in the Hindu Kush 
crustal and mantle zones warrants two separate investigations (Wheeler 
et al., 2005). This study is aimed at understanding the Hindu Kush 
seismicity by comparing the shallow (0-70 km) and intermediate depth 
(70-300 km) earthquakes using various earthquake histograms, time- 
series plots, b-values, depth-dependent b-values, low and high b-value 
zones and z-values. 

2. Seismotectonics of the Hindu Kush 

The present work uses shallow (0-70 km) and intermediate 
(70-300 km) depth earthquakes which occurred in the area bounded by 
34.35-39°N and 68.4-74.47°E (Fig. 2). The magnitude interval con- 
sidered was 4.08-7.7 during the period of 1900-2007 (Rehman et al., 
2014). Wheeler et al. (2005) discussed the relationship of the faults to 
the major tectonic zones of the Hindu Kush. The location of faults is 
depicted in Fig. 2 as black lines. The Darvaz fault is seismically active, 
lying in Afghanistan and Tajikistan. Trifonov (1978) suggested that the 
Darvaz fault is a region of high seismic slip toward the Afghanistan side, 



Fig. 3. Time histograms of (a) shallow earthquakes (zone A; 0-70 km) 
and (b) intermediate earthquakes (zone B; 70-250 km). 
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Fig. 4. Histograms of the magnitudes of earthquakes with (a) zone A 
and (b) zone B. 


and earthquakes concentrate along this part of the fault. The other 
active faults are the Central Badakhshan fault and the Konar fault. The 
central part of the Central Badakhshan fault contains intense seismic 
activity extending to depths up to 240 km. Earthquakes of magnitude 
7.0 and greater are primarily concentrated in this part of the fault 
(Table 2). Fig. 2 shows that intermediate depth earthquakes have S- 
shaped seismicity, while shallow depth earthquakes are primarily 
concentrated in the North Afghan Platform (69°N, 36.4°E), the junction 
of the Konar fault and MMT, Pamir (73°N, 37°E) and in the region near 
the Central Badakhshan fault. 

3. Methods 

Two important seismological parameters used in frequency-magni- 
tude relation are the b-value and magnitude (M) of an earthquake. The 
b and M of the Gutenberg and Richter (1944)-GR relationship is re- 
presented by 

LogN = a-bM (1) 

A number of b-values (the slope of the GR relation) studies have 
been done in previous years (e.g., Mogie, 1962; Wyss, 1973; Seeber 
et al., 1980; Urbancic et al., 1992; Schorlemmer et al., 2005; Khan and 
Chakraborty, 2007; Sarkarinejad et al., 2013; El-Isa and Eaton, 2014; 
Maden and Ozturk, 2015; Singh and Singh, 2015). These findings 
suggest that the value of b may be related to tectonic, material 


heterogeneity and stress level. The higher b-value is due to weak and 
heterogeneous rocks, lower stress and numerous smaller earthquakes, 
while the smaller values of b indicate an increase in stress and homo- 
geneity. 

The robust method of maximum curvature is used to evaluate the 
temporal variation in the magnitude of completeness, and its un- 
certainty from the GR relation (Wiemer and Wyss, 2000; Woessner and 
Wiemer, 2005; Fig. 9). The method of maximum curvature suggests that 
the completeness magnitude (Me) of Mw 5.0 seems reasonable for the 
region of Hindu Kush. Zarifi and Havskov (2003) used a threshold value 
of 4.8 Mb (body wave magnitude) for the Hindu Kush. This value of Mb 
is approximately equal to Mw 5.0 using the empirical global magnitude 
conversion relations of Scordilis (2006). The earthquake catalogue for 
the Hindu Kush used is drawn from the magnitude homogenized 
earthquake catalogue of Rehman et al. (2014). The homogeneous 
earthquake catalogue is the primary input in the study of seismicity 
variation (Habermann, 1987; Burton, 1990; Katsumata, 2011). Rehman 
et al. (2014) purified various existing earthquake catalogues (e.g., In- 
ternational Seismological Centre (ISC), EHB (Engdahl-Hilst-Buland; 
Engdahl et al., 1998), U.S. Geological Survey’s National Earthquake 
Information Centre (NEIC), Global Seismic Hazard Assessment Pro- 
gramme (GSHAP; Zhang et al., 1999), ABE - Abe and Kanamori (1979) 
and SSR - Kondorskaya and Shebalin (1982)), and used appropriate 
magnitude conversion relations (Scordilis, 2006) to obtain a homo- 
genized earthquake catalogue for Pakistan. All magnitude values 
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Fig. 5. Histograms in (a) and (b) show depth difference with respect 
to the number of events for A and B seismicity, respectively. 
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Fig. 6. Plots of earthquake magnitude versus time for (a) 
zone A and (b) zone B earthquakes. 
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Fig. 7. Depth variation of the seismicity for (a) zone A and 
(b) zone B. 



discussed in the current study are expressed in moment magnitude 
(Mw) scale using a homogenized earthquake catalogue of Pakistan after 
Rehman et al. (2014). The homogenized catalogue for Pakistan has 
been produced using the abovementioned catalogues data taken from 
the British Geological Survey (BGS; The BGS World Seismicity 
Database, 2004) and National Engineering Services Pakistan (Pvt) 
Limited (NESPAK; NESPAK, 2006). Table 1 summarizes the basic sta- 
tistics of the implemented catalogues for the Hindu Kush. 

When analyzing the spatial-temporal behavior of earthquakes, it is 
most important to consider a homogeneous earthquake catalogue 
(Burton, 1990; Burton et al., 2004; Bayliss and Burton, 2007; Yadav 
et al., 2011; Das et al., 2012; Rehman et al., 2014). However, there is 
often a need to describe the parameters of the earthquake catalogue in 
error terms. In this regard, Maggi et al. (2000), Waldhauser and Schaff 
(2007) and Camanni et al. (2016) developed useful methods to enable 
errors in regional and global catalogues. Maggi et al. (2000) used the 
method of synthetic seismogram to study the depth distribution of 
earthquakes. Waldhauser and Schaff (2007) developed a double-dif- 
ference method to reveal improvements in global event locations. 
Camanni et al. (2016) used a statistical method of collapsing for stan- 
dard errors in hypocenter locations. The implemented earthquake cat- 
alogue for the Hindu Kush is dominated by major available sources of 
ISC and EHB data (Table 1). ISC routinely update the global earthquake 
parameters data and provide the most comprehensive available earth- 
quake catalogue for seismic hazard studies (Burton et al., 2004; 
Waldhauser and Schaff, 2007; Yadav et al., 2011). ISC reanalyze and 
expand the existing NEIC catalogues using additional seismic phase 
data. The updated catalogue of EHB (Engdahl et al., 1998; akl35 
model) is generally considered to be a modern and accurate catalogue 


that is available for global seismic studies due to the significantly im- 
proved hypocenter determination (Bondar et al., 2004; Pesicek et al., 
2010). Elements that serve to control the EHB model are subsequently 
identified. For example, relocation of well-constrained hypocenter lo- 
cations, phase arrivals of P, S and PKP, teleseismic depth phases pP, 
pwP and sP, Ellipticity, station corrections and weighting pattern using 
phase variance as a function of distance. The approximate depth error 
in the EHB catalogue is 10 km, which is slightly lower than the ISC 
catalogue (Pesicek et al., 2010). 

The present analysis uses the method of maximum-likelihood as the 
most appropriate and more reliable approach for estimation of the b- 
value than the least square technique (Aki, 1965; Utsu, 1965; Bender, 
1983; Woessner and Wiemer, 2005; Singh and Singh, 2015). b- value 
analysis using the method of maximum likelihood is given by the fol- 
lowing relation: 

b = log 10 e(l/M-M min ) (2) 

where M and M min are the average magnitude and threshold magni- 
tude, respectively. 

Z-statistics (Wiemer and Wyss, 1994; Kawamura et al., 2014) for 
comparing the spatio-temporal seismicity rate changes of the areas is 
given below: 

z(0 = , 

ni + o) 

where Xi, oi and ni are the mean seismicity rates, the variances and the 
number of earthquakes in two different time periods, respectively. 
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Fig. 8. Cumulative seismicity curves for (a) zone A and (b) 
zone B. 
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4. Results and discussion 

In this study, two zones, A (shallow earthquakes) and B (inter- 
mediate earthquakes) are compared and analyzed using (1) histograms 
of earthquakes (frequency of earthquakes vs time, magnitude and 
depth), (2) time series plots, (3) the b-values (depth-dependent-b and 
the resulting high b-values zones), and (4) z values, which are shown in 
Figs. 3-5, Figs. 6 and 7, Figs. 8-10 and Figs. 11 and 12, respectively. 

Figs. 3-5 show the relevant time, magnitude and depth histograms 
of the seismicity data. From the time histograms of shallow (Fig. 3a; 
0-70 km) and intermediate (Fig. 3b; 70-250 km) depth earthquakes in 
the Hindu Kush, it is apparent that seismicity data are divided into two 
main types of time periods (pre- and post- 1960). This is simply due to 
the large number of earthquake events that were reported in the post- 
1960 period. Both histograms show a dramatic rise in the number of 
earthquakes reported during the post- 1960 period. This increase is due 
to the international seismograph networks of WWSSN with an im- 
proved and developed earthquake monitoring system during the 1960s 
(Burton et al., 2004). Fig. 3a and b compares the pre- and post-1960 
distribution of shallow and intermediate depth earthquakes in the 
Hindu Kush. It helps to explain that the increase in the zone B seismicity 
is much higher than the zone A. Zone B shows a higher frequency end 
(approximately 240) than that of A (approximately 60). Similarly, the 


frequency of zone B events is higher (between 0 and 25) than the zone A 
(between 0 and 8) pre-1960. Comparison of zone A and B events is 
further evident from the post- 1960 data. Confirmation for this finding is 
provided by the significant increases in the reported earthquakes (more 
than 200) around two prominent peaks of 1996 and 2001 (Fig. 3b) 
compared to the 1983 and 1988 peaks in Fig. 3a (70 or greater). 
Magnitude histograms of zones A and B are shown in Fig. 4a and b, 
respectively. Both histograms divide the earthquake data into two 
prominent magnitude ranges: (a) 4.08-5.1 in Fig. 4a, b and (b) 5. 2-7.0 
in Fig. 4a and 5. 2-7.7 in Fig. 4b. On the far right of the histograms, we 
see low bars until the magnitude 7.0 in zone A, while zone B seismicity 
displays earthquakes beyond magnitude 7.5. Only one maximum peak, 
at magnitude 4.7, is identified in zone A, where earthquake frequency 
has reached to 100. Conversely zone B has three prominent peaks at 
magnitudes 4.4, 4.7 and 5.0, where earthquakes are more than 270. The 
frequency of earthquakes decreases beyond magnitudes 5.1 and 5.4 in 
zone A and zone B, respectively. Zones A and B depth histograms are 
presented in Fig. 5a and b, respectively. In addition to one prominent 
peak that has high frequency ( > 250 earthquakes) at a depth of 33 km 
(Fig. 5a), which is the ISC default depth for unknown earthquakes depth 
(Bolton et al., 2006), there are trends of lower frequency earthquakes in 
a wider range of depth. Three other peaks of earthquake activity are 
also observed at 1, 10 and 15 km, where a frequency of more than 30 
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Fig. 9. Cumulative Frequency - Magnitude Distribution of (a) zone A using 1964-2007, 
1900-2007, (e) zone B using 1900-2007, (f) whole of the Hindu Kush using 1900-2007. 


seems to occur. All other depth bars display the frequency of earth- 
quakes below 20. Frequency of earthquakes increases in the zone B 
histograms, which show considerable variation in the geographical area 
of the Hindu Kush (Fig. 5b). Furthermore, it is observed in Fig. 5b that 
the two depth bar ranges of 88-120 km and 200-220 km have a higher 
frequency than the remaining bars. 

The magnitude distribution in zones A and B of earthquakes as a 
function of time is shown in Fig. 6a and b. Both figures show that the 
earthquakes magnitude reported in the Hindu Kush region are primarily 
after 1963. Most of the earthquake magnitudes were reported in the 
time span of 1964-2007. Reported magnitude events greater than 6.5 
are likely complete since 1900. Similarly, earthquakes as a function of 
depth increases during the same time period of post- 1963, and the time- 
depth plots of zones A and B seismicity data are displayed in 
Fig. 7a and b, respectively. Furthermore, similarities and dis-similarities 
can be seen during pre- and post- 1963 in both plots. 

Fig. 8a and b shows the cumulative number of earthquakes vs time 
in zones A and B earthquakes in the Hindu Kush, respectively. Sig- 
nificant changes are clearly observed in the cumulative number of 
earthquakes as a function of time after 1963. Approximately 70 and 250 
earthquakes are reported at the starting date of 1964 in A and B seis- 
micity curve, respectively. Similarly, the maximum reported earth- 
quakes are approximately 820 and 3020 in zones A and B earthquake 



Magnitude 




(b) zone B 1964-2007, (c) whole of the Hindu Kush using 1964-2007, (d) zone A using 


plots, respectively. Time histograms (Fig. 3a and b), time variation of 
seismicity (Fig. 6a and b) and cumulative seismicity curves 
(Fig. 8a and b) provide the same and consistent time interval of 
1964-2007 for the study of b-value calculation as shown in 
Fig. 9a and b. The b-value, its standard deviation and a-value of the 
Gutenberg-Richter relation are calculated using the maximum like- 
lihood method because it is a more robust estimation than least-square 
(Aki, 1965). The time interval of 1964-2007 is considered for b-values 
calculation of zones A and B. Fig. 9a and b shows the plots of cumu- 
lative number of earthquakes against magnitudes of zones A and B 
earthquakes in the Hindu Kush, respectively. The b-value is calculated 
as 1.32 ± 0.07, and Me is calculated as 5 for shallow seismicity 
(Fig. 9a). For intermediate depth earthquakes, shown in Fig. 9b, the li- 
vable is calculated as 1.07 ± 0.03 using the same completeness value 
of Me = 5 mentioned earlier for Fig. 9a. The b-value calculation sug- 
gests that the zone A seismicity (b-value = 1.32) is clearly separated 
from the zone B earthquakes (b-value = 1.07). The finding indicates 
that large magnitude earthquakes occurred in the area of high stress of 
low value (intermediate depth), which correlates well with the current 
scenario (Figs. 4 and 6). The b-value for whole Hindu Kush is found to 
be 1.06, which is shown in Fig. 9c. This value is close to zone B, which 
is due to the greater portion of seismic activity concentrated in the 
intermediate depth. Calculation of the b-value for zone A, B and the 
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b-value 


whole Hindu Kush is made using the time interval of 1900-2007. B- 
values for zone A, zone B and the whole Hindu Kush are shown in 
Fig. 9d, e and f, respectively. Again, the b-value for zone B (0.802) is 
lower than the b-value for zone A (0.907) and close to the value for the 
whole Hindu Kush (0.852). The results also reflect the consistent high 
seismic area of intermediate earthquakes. The b-value range of 0.5-1. 5 
indicate the tectonic earthquakes environment and the general value is 
approximately 1.0 (Reasenberg and Jones, 1989). Gutenberg and 
Richter (1949) estimated the global range of b-value from 0.45 to 1.5 
for shallow, intermediate and deep earthquakes. Miyamura (1962) 
suggested a b-value range of 0. 7-1.0 for the region of the Circum-Pa- 
cific and Alpine orogenic belts. Zarifi and Havskov (2003) determined 
the value of b-value in the Hindu Kush is 1.43 using intermediate depth 


earthquakes. 

Fig. 10 represents the average b-value with standard deviations 
(horizontal bars) versus the depth of earthquakes in the Hindu Kush. 
Fig. 10a clearly identifies the ISC default depth of 33 km, where 
earthquake depth is not available (Bolton et al., 2006). The b- values 
show a high peak around the depth of 33 km. The b-value range of 
1-1.2 is found above the default depth while it increases from 1.2 to 1.5 
below 33 km. b-values smaller than 1 are identified in zone B events 
plot (Fig. 10b). Considerable variation can be seen in the b-value vs 
depth graph. According to Singh and Singh (2015), b less than 0.9 
describes an area of low values for the Hindu Kush and b greater than 
0.9 are suggested as high. Following this finding, three zones are as- 
sociated with the occurrence of intermediate-depth earthquakes having 
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Fig. 11 . Low b-value areas of intermediate earthquakes in the Hindu 
Kush. 


b-values lower than 0.9. Zone 1 contains the depth range of 
110-160 km. Zone 2 provides low b-values approximately 175 km, and 
zone 3 considers depth variation ranging from 205 to 240 km. These 
show considerable variation in the geographical region of the Hindu 
Kush, as shown in Fig. 11. The clear distinction occurs in the three low 
b-values regions indicated by blue, red and green colors. More variation 
occurs in the seismic zone with a depth range of 110-160 km, which 
extends to the northeast of the Hindu Kush. Clusters of earthquakes 
occurred at 71-72°E and 36-37°N and contrasts with northeastward 
and southwestern earthquake events. This zone generated earthquakes 
with magnitudes up to 7.0. Few earthquakes of low b-values at a depth 
of 175 km (red color) are shown near the above clustering group of 
earthquakes. This zone of few earthquakes is characterized by a mag- 
nitude of 5.3 or less. Zone 3 earthquakes (green color) are located at 
71-71.5°E and 36.25-36. 73°N, and the zone become narrow towards 
southwest. Large earthquake magnitudes of 7. 0-7. 7 occur in this zone 
of low b-value. Chatelain et al. (1980) and Zarifi and Havskov (2003) 
provided the general interpretation of intermediate depth earthquakes 
of the Hindu Kush and suggested that the earthquakes occur in sepa- 
rated clusters with aseismicity gaps between the clusters. 

Fig. 12a and b shows the z-maps for zones A and B earthquakes, 
respectively. The z-map represents the spatial distribution of the dif- 
ferences of the seismicity changes and it compares the seismicity rates 
from 1900 to 1965 to those from 1965 to 2007. The blue color (de- 
creased Z-value) and red color (increased Z-value) represent increases 
and decreases in the seismicity rate changes, respectively. The area 
shown in Fig. 12a is characterized by low seismicity rate changes except 
the default depth region approximately 36.5°N and 71.5°E. The Z-value 
map for zone B earthquakes (Fig. 12b) shows that there is consistent 
and significant increase in the seismicity rate changes in the three high 
seismic activity zones as described in Fig. 11. Furthermore, the region 
of high seismicity rate change has experienced significant earthquakes 


in zone 3 with magnitudes of 6.6 (1911, 1921 and 1949) and 7.7 (1940 
and 1949) as shown with the yellow stars (Table 2). 

5. Conclusions 

In this study, two separate investigations were conducted to com- 
pare and analyze the distribution of shallow (0-70 km) and inter- 
mediate (70-300 km) depth events in the Hindu Kush using various 
histograms, time-series plots, b and z values. The time variation of 
seismicity (Fig. 3a and b) shows that more earthquakes are reported 
during 1964-2007. This shows high seismicity periods after 1980 in 
case of shallow earthquakes and after 1964 in case of intermediate 
depth earthquakes. The relation between magnitude and number of 
earthquakes has been plotted in a magnitude histogram, which re- 
presents high percentage of magnitude reported for intermediate depth 
earthquakes than shallow seismicity. Furthermore, intermediate depth 
seismicity is characterized by high magnitude earthquakes (greater 
than 7.5) than shallow earthquakes (only one earthquake with magni- 
tude 7.0). The depth histograms indicate that the frequency of earth- 
quakes is dependent on depth. Approximately 79% of all earthquakes 
with a depth greater than 70 km are observed, while shallow seismicity 
region contains 21% of earthquakes. The frequency increases con- 
siderably when the depth is below 70 km. From the histograms and 
time series plots, we conclude that the percentage of reported earth- 
quakes is increased during the post-1964 WWSSN period. 

The b-value for shallow and intermediate earthquakes is calculated 
during 1964-2007. Shallow seismicity region represent relatively 
higher value of b and shows lower seismicity region. The region of 
intermediate depth earthquakes has lower b-value of 1.07 and hence 
shows higher seismicity. This value of 1.07 is close to b-value for the 
entire Hindu Kush. The same seismicity pattern is determined with a 
catalogue period of 1900-2007. Seismic b-value within shallow 
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Fig. 12. Seismicity changes for (a) zone A and (b) zone B. 


earthquake zones is 0.907. b-values of 0.852 and 0.802 are observed for 
the intermediate earthquake zone and of the entire Hindu Kush, re- 
spectively. Furthermore, the intermediate earthquake region exhibited 
a variation in b- value and the various zones of b- value less than 0.9 are 
shown and plotted. Three lower b-values (less than 0.9) areas are 
suggested: (1) depth range of 110-160 km, (2) depth of 175 km, and (3) 
depth variation ranged from 205 to 240 km. These results and the 
earthquake distribution are in agreement with the z-models. There is a 
consistent increase in the seismicity rate changes in the suggested low 
b-value areas. 

In this study, we show the importance of comparing and separating 
seismicity of the Hindu Kush into two regions of shallow and inter- 
mediate depth using an attractive and advanced seismological tool, 


ZMAP. It is shown that the depth range of 70-300 km is more seismi- 
cally active than the shallow depths. Today, many algorithms are 
available in the study of seismicity; however, emphasis should be 
placed on the need to consider the associated errors within the im- 
plemented earthquake catalogue. 
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Table 2 

Earthquakes of magnitude 7.0 and greater in the Hindu Kush. 


Year 

Month 

Day 

Latitude 

Longitude 

Depth (km) 

Mw 

1902 

10 

6 

36.50 

70.50 

200 

7.3 

1903 

4 

19 

37.00 

71.00 

160 

7.0 

1907 

12 

25 

36.50 

70.50 

240 

7.0 

1908 

10 

24 

36.50 

70.50 

220 

7.1 

1909 

7 

7 

36.50 

70.50 

230 

7.2 

1911 

7 

4 

36.00 

70.50 

190 

7.6 

1917 

4 

21 

37.00 

70.50 

220 

7.3 

1921 

5 

20 

36.00 

70.50 

220 

7.2 

1921 

11 

15 

36.50 

70.50 

215 

7.6 

1922 

12 

6 

36.50 

70.50 

230 

7.2 

1929 

2 

1 

36.50 

70.50 

220 

7.4 

1933 

1 

9 

36.40 

69.60 

180 

7.1 

1937 

11 

14 

36.50 

70.50 

240 

7.2 

1939 

11 

21 

36.30 

70.60 

220 

7.3 

1940 

3 

4 

36.60 

70.60 

99 

7.7 

1943 

2 

28 

36.50 

70.50 

210 

7.1 

1944 

9 

27 

39.00 

73.50 

40 

7.0 

1949 

3 

4 

36.00 

70.50 

230 

7.7 

1956 

11 

14 

36.70 

70.80 

190 

7.0 

1960 

2 

19 

36.80 

71.00 

190 

7.1 

1962 

7 

6 

36.80 

70.10 

210 

7.4 

1965 

3 

14 

36.40 

70.72 

210 

7.6 

1966 

6 

6 

36.42 

71.14 

216 

7.2 

1971 

8 

4 

36.39 

70.74 

202 

7.2 

1974 

7 

30 

36.37 

70.74 

211 

7.3 

1983 

12 

30 

36.39 

70.71 

214 

7.4 

1985 

7 

29 

36.17 

70.87 

100 

7.4 

1993 

8 

9 

36.33 

70.88 

211 

7.0 

2002 

3 

3 

36.43 

70.44 

209 

7.4 
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